!   _______________________________________________________________________________
!    _______/\\\\\\\\\_______/\\\\\\\\\______/\\\\\\\\\\\\\\__/\\\\\\\\\\\\\\\______
!     _____/\\\////////_____/\\\////////_____\/\\\//////////__\/\/////////////________
!      ___/\\\/____________/\\\/______________\/\\\____________\/\\\________/\\\\_______
!       __/\\\____________/\\\_________________\/\\\\\\\\\\\\\__\/\\\_______\////________
!        _\/\\\___________\/\\\_________________\/\\\/////////___\/\\\________/\\\\_______
!         _\//\\\__________\//\\\________________\/\\\____________\/\\\_______\////________
!          __\///\\\_________\///\\\______________\/\\\____________\/\\\______/\\\\_________
!           ____\////\\\\\\\\\____\////\\\\\\\\\___\/\\\____________\/\\\\\\\\\////__________
!            _______\/////////_______\/////////_____\///_____________\/////////_______________
!            __________________________________________________________________________________
!----------------------------------------------------------------------------------------------
!>  subroutine input_files_cgns
!>  last edit 2016-05-24
!>  last edit by liuxz
!----------------------------------------------------------------------------------------------
    subroutine initial
!=======================================================================
        use global_parameter
        !>
        !>
		implicit none
        real(kind= dprec) t1,t2,u1
        !>
        !>
        !>
        if(mach .lt. 0.0 )then
            write(*,'(" ccfdv3.0 :: the mach number you input is negative and the mach= ",e12.6)') mach
            xmach = abs(mach)
        else if(mach .eq. 0.0)then
            write(*,'("ccfdv3.0:: the mach number you input is 0.0")')
            xmach = 0.6
        else if (mach .gt. 1.0)then
            write(*,'("ccfdv3.0:: the flow is supersonic,")')
            write(*,'("and the all-speed computing capacity have not fully develop yet!!!!")' )
            write(*,'("the furture the ccfd3.0 will have the all-speed computing capacity")')
            xmach = mach
        else
            xmach = mach
        end if

        !>
        !>
        rmre = xmach/reue
        !>
        r0   = 1.
        p0   = r0/gamma
        u0   = xmach*cos(alpha)*cos(beta)
        w0   = xmach*sin(alpha)*cos(beta)
        v0   = -xmach*sin(beta)
        !>
        !>
        u1 = 0.5e0*(u0**2+v0**2+w0**2)
        t1  = p0/r0/(gamma-1.0)
        t2  = r0*(t1+u1)
        h0  = (t2 + p0)/r0

        return
    end subroutine initial
